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Hawking Radiation and Non-equilibrium Quantum Critical Current Noise 
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The dynamical scaling of quantum critical systems in thermal equilibrium may be inherited in the 
driven steady-state, leading to universal out-of-equilibrium behaviour. This attractive notion has 
been demonstrated in just a few cases. We demonstrate how holography - a mapping between the 
quantum critical system and a gravity dual - provides an illuminating perspective and new results. 
Non-trivial out-of-equilibrium universality is particularly apparent in current noise, which is dual to 
Hawking radiation in the gravitational system. We calculate this in a 2-dimensional system driven 
by a strong in-plane electric field and deduce a universal scaling function interpolating between 
previously established equilibrium and far-from-equilibrium current noise. Since this applies at all 
fields, out-of-equilibrium experiments no longer require very high fields for comparison with theory. 

PACS numbers: 



Relatively few general principles are known that govern 
the behaviour of quantum systems driven out of equilib- 
rium [IHl] ■ Those discovered do not constrain the system 
as tightly as those of equilibrium thermodynamics. For 
quantum critical systems, there may be additional prin- 
ciples leading to universal out-of-equilibrium behaviour. 
Quantum critical systems display dynamical scaling in 
thermal equilibrium [5]. Since non-equilibrium steady 
states are largely determined by dynamics, one might an- 
ticipate that this scaling is inherited by the driven steady 
state. This has been seen in a few cases [BUTT]. However, 
analysis of out-of-equilibrium behaviour is fraught with 
conceptual and calculational difficulties. New insights 
and approaches are required. 

Holography provides one such approach [TH [T3]. It 
utilizes a mapping between gauge theory in d dimensions 
and string theory in d -|- 1 dimensions; the classical grav- 
ity saddlepoint of the string theory corresponding to a 
strongly-correlated quantum state of the gauge theory. 
In this way, holography encodes quantum many-body 
correlations in a classical, gravitational metric, allowing 
complementary physical insights and calculational tools 
to be brought to bear [T^HTH]. 

Here, we apply the holographic technique to the study 
of non-linear current noise - dual to Hawking radiation - 
at a two-dimensional z = \ quantum critical point driven 
by an in-plane electric field. We find a noise power 
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Sj ^ AaT^ with nksT^ = [(tt/cbT)"' 

(1) 
interpolating between Johnson- Nyquist noise at low fields 

and previous results at high-fields [5]. Whilst driven 
out-of-equilibrium steady-state distributions are gener- 
ally very different from those in equilibrium, holography 
reveals that the out-of-equilibrium fluctuations of a quan- 
tum critical system may be precisely thermal. 
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FIG. 1: Sketch of the AdS geometry with boundary current 
fluctuations: The AdS brane geometry is modified by the ap- 
plication of an m-plane electric field at the boundary. Hawking 
radiation from the modified horizon propagates to the bound- 
ary where it induces current fluctuations. These fluctuations 
satisfy an equilibrium-like fluctuation dissipation relation at 
a temperature T, — \T*' -\-h}c*'e'^E^ /{■KkB)'^Y , interpolating 
between established low- and high-field limits. 



Our calculation is carried out in a probe-brane limit 
with boundary conditions corresponding to the applica- 
tion of the in-plane electric field [IT . This field induces a 
modified metric [.19j on the probe brane with a horizon at 
temperature T* [38]. We show explicitly that Hawking 
radiation from the emergent horizon generates current 
fluctuations in the boundary that are thermal at tem- 
perature T* (see Fig. 1), a conclusion consistent with 
general gravitational principles and with potentially pro- 
found implications for out-of-equilibrium quantum criti- 
cality 

These results are compared to those for the bosonic 
superconductor to insulator transition, embodied in the 
Bose-Hubbard model. This is amongst the quantum crit- 
ical systems that we have under the best analytical con- 



trol [201 121] ! making it ideal to test new ideas and tech- 
niques. Early work on out-of-equilibrium quantum crit- 
icality considered this system |71 [S]. Pioneering holo- 
graphic studies have reproduced dualities of the equilib- 
rium magneto-thermo-electric response [THl [22] and the 
steady-state out-of-equilibrium current [H [11] . These re- 
sults are a useful benchmark and demonstrate the utility 
of the holographic approach. 

We begin by outlining how the non-equilibrium steady- 
state response and current noise are obtained in a 1/A^- 
expansion of the Bose-Hubbard model. We follow with a 
heuristic discussion of the holographic mapping and how 
it can be applied to out-of-equilibrium steady-states us- 
ing the probe brane limit. We pay particular attention 
to the similarities between the trick used in the 1/A^ ex- 
pansion and the probe brane limit. Next, we detail the 
calculation on the gravity side and end with a critical 
discussion of the results and their implications. 

Out- of- Equilibrium Noise in the Bose-Hubbard Model: 
The continuum field theory describing the two- 
dimensional Bose-Hubbard model is essentially a Klein- 
Gordon model with a (f>* interaction, minimally coupled 
to an in-plane electric field [201 . The Klein-Gordon field 
supports normal modes of positive and negative charge. 
Charge transport in this model can be described using a 
Boltzmann equation for the population of positively and 
negatively charged normal modes [21| , controlling the in- 
teraction within an e-expansion. The presence of both 
charges allows scattering to relax current even though 
it conserves both energy and momentum. Formally, the 
model has infinite thermal conductivity, k, at leading or- 
der in e and finite electrical conductivity, a. 

As the electric field is increased and the system 
driven away from linear response, new processes be- 
come important and account must be taken of heat 
flow. The strong electric held produces pairs of positive 
and negative modes from the vacuum by the Schwinger 
mechanism/Landau-Zener tunnelling j7]. At the same 
time. Joule heating pumps energy into the system. In 
order to achieve a steady state, there must be a heat 
sink. One way is to consider a finite system so that ex- 
cess heat and particles can be transported to the bound- 
aries. As emphasized by Tremblay [23], heat transport 
requires gradients in free-energy density. If this vari- 
ation is to be small, and the response independent of 
the size and geometry of the sample, there is a limit on 
size. However, for particle-hole symmetric systems, the 
Wiedemann-Franz law does not hold and thermal con- 
ductivity can be much larger (in appropriate units) than 
electrical conductivity. In this case, a very small ther- 
mal gradient is required to support the flow of heat to 
the edge of the sample and maintain a steady state. In- 
deed, for the Bose-Hubbard model, the leading thermal 
conductivity is infinite and the system size may be taken 
to infinity provided that we take k to infinity first |22j . 
Both the 1/N expansion and the probe brane calcula- 
tion take advantage of these arguments for the existence 
of a steady-state out-of-equilibrium distribution and cal- 



culate the out-of-equilibrium distribution directly for an 
infinite system. 

In the 1 /TV-expansion [8 , the electric field is coupled to 
just one of the components of an 0{N) vector field. The 
A^— 1 components that are not coupled to the electric field 
effectively act as a critical bath for the field-coupled com- 
ponent. This allows the steady-state out-of-equilibrium 
conductivity to be calculated in a spatially uniform, infi- 
nite system. A Boltzmann equation is developed for the 
field-coupled component of the 0{N) field. Integrating 
over momentum, one recovers a Boltzmann equation for 
the total current; 



dt 
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where the first term on the right hand side has its origin 
in Schwinger pair production and the second term de- 
scribes the decay of current carried by the field-coupled 
component to the remaining A^ — 1 components. The con- 
stants a and b characterise these two processes [33] • Inte- 
grating this equation, one obtains a steady-state current 
j — s^Ed/2 giving a conductivity cr = ^ in two dimen- 
sions. Note that in this 1/N treatment, the steady-state 
distribution is not simply thermal, a point that we shall 
return to later. 

The current noise is determined by developing a 
Boltzmann-Langevin [231 di] description of fluctuations 
in the distribution function. The underlying physical in- 
sight is that the scattering processes may be treated as 
independent and so obey Poisson statistics. Fluctuations 
in the distribution function may then be approximated by 
adding Gaussian noise to the Boltzmann equation with 
variance proportional to the mean scattering rate. Inte- 
grating this Boltzmann-Langevin equation results in the 
following equation for the current fluctuations: 
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Integrating this equation and considering the behaviour 
on long timescales, one obtains the large- -E (zero temper- 
ature) limit of the non-equilibrium current noise given by 
Eq.([l]). 

Heuristic Picture of the AdS Realization: As discussed 
above, holography uses a mapping between a critical 
gauge theory in d dimensions and a string theory in d-l- 1. 
The gauge theory may be thought of as describing or liv- 
ing on the boundary of the gravitational system with 
the extra dimension playing the role of a renormaliza- 
tion group scale. The classical gravity saddlepoint of 
the string theory encodes a strongly-correlated quantum 
state of the gauge theory. In particular, the interactions 
and scaling are embodied in the metric. The horizon in 
the presence of a black hole has particular importance in 
what follows. As usual, its area determines the entropy of 
the black hole, whose temperature is given by the surface 
gravity. This temperature has the same interpretation in 
both the gauge theory and gravity. 



We model non- linear transport using the D3n/D5m 
brane intersection [26l [27] : we study the dynamics of M 
D5 branes in the background of ND3 branes in the probe 
brane limit N 3> M. The specific configuration is sum- 
marized in table U The D3 branes lead to the formation 
of an horizon in the background metric modelling the ef- 
fect of a thermal bath. Non-linear transport is studied 
through the dynamics of gauge fields on the D5 brane 
with boundary conditions corresponding to the applica- 
tion of an in-plane electric field. This permits us to study 



TABLE I: The D3n/D5m brane intersection. We la- 
bel directions occupied by the branes by • and directions in 
which they are localised by o. The field-theory directions are 
a;", x^,x^ = t, X, y. In addition the D5 wraps the RG direction 
^ — u as well as an 5"^ ^ '^^ 

,6 



X 



C S^, parametrised by two angles 
x",x" = til, (pi. This system realises an out-of-equilibrium 
steady state in 2 + 1 dimensions. 
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out-of-equilibrium steady states in a very similar way to 
the 1/iV trick used in the condensed matter calculations 
of current noise. Since N ^ M, gauge fields on the D5 
brane do not modify the metric. Just as the single com- 
ponent of the 0{N) field coupled to the electric field does 
not modify the distribution of the TV — 1 modes corre- 
sponding to the critical bath from the condensed matter 
point of view. 

The application of the electric field on the boundary 
leads to a crucial modification of the metric induced on 
the D5 brane. A new horizon appears in this reduced 
metric with an area that corresponds to the effective tem- 
perature T* in Eq.Q. A major result of this is that the 
conductivity and current fiuctuations are related in a way 
that mimics the equilibrium fiuctuation-dissipation rela- 
tion at the effective temperature T*. Such a relation ap- 
pears to be inevitable in the gravity dual. As we discuss 
below, the consequences of this are potentially profound. 

Details of the AdS Calculation: Our calculation, in- 
volves solving the equations of motion for gauge fields on 
the D5 brane. The steady-state solution corresponds to 
a previous calculation that allows determination of the 
out-of-equilibrium conductivity |Tl]. We explicitly prop- 
agate thermal fluctuations from the horizon on the D5 
brane to its boundary and derive a Langevin equation 
for gauge fluctuations and hence current fluctuations on 
the boundary. The dynamics of the D3n/D5m brane 
intersection is captured by the action 



Sd, = -A^s f (fi^-det (5i„d + F) + ^wz , 



(2) 



is the gauge fleld strength on the brane in string units. 
We do not specify the Wess-Zumino term 5'wz further 
as it does not contribute to the configurations that we 
study below. The normalisation A/Js = MT^ where T5 is 
the usual Z35-brane tension. We study this system in the 
background of the non-extremal DZ brane metric, which 
we take to be 
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■'f{u)dt^ + dx^) 



---^ + R^d^l 
f{u)u^ 

(3) 
with f{u) = 1 — ^ , and dilg is the metric on a round five 
sphere. By studying the regularity of the Euclidean ge- 
ometry at Uh or otherwise, one sees that the background 
geometry (the bath) has a temperature T = Uh/i'n'R'^). 

Applying an external electric field in the x direction is 
accomplished by choosing the gauge potential 



(4) 



A= {-Et-A^{u)) dx, 



where £^ is a dimensionless electric field. The geometry of 
the brane is obtained by solving the equations resulting 
from the action ^. We choose the static gauge on the 
brane. The degrees of freedom are a single angle 9(u) de- 
scribing the position of the S*^ on the S^ and the radially 
dependent function Ax{u). An important simplification 
allows us to set 9{u) = 0, corresponding to the case of 
coupling massless charged degrees of freedom to the field. 
The field Ax{u) is a cyclic coordinate and hence has 
first integral 
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where C is an integration constant related to the current 
G = j /Nz- Requiring that the expression be real for all 
values of u demands |28j that denominator and numera- 
tor cross zero at the same point u := m*. This determines 
M^ = u^ -f E^R'^. The special point u* on the brane co- 
incides with the horizon of the background black hole iff 
E = 0. This matching procedure also relates the current 
to the applied field, yielding (see [TT] ) 
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M-.E 



C^(2+l) 



E, 



(6) 



where gind is the metric induced on the D5 brane, and F 



i.e. non-linear conductivity in d = 3 is a constant IHJ'T]- 
Brane Horizon, Effective Temperature: This result, 
and its higher-dimensional generalisations, can further 
be interpreted in terms of an effective geometry on the 
brane, seen by the fluctuations of the gauge field. Let us 
perturb A — > A + a, where 

a ^ a"{u,t)dx + a {u,t)dy + a^{u,t)dt . (7) 

Then the quadratic fiuctuation action becomes 



5^2) = JVr, f dudtV^a''^ (daJdbJ + Z{u)^daJdbJ + Z*{u)daa'dba* 



(8) 



where A/Js results from integrating out the angular di- 
rections, and for the case of no spatial dependence also 
the x,y directions. Z^{u) are conformal rescaling factors 
whose exact expressions wc do not need. Upon plug- 
ging in the conductivity ^ we find that Z-^{u) = 1. By 
choosing a new time variable r we can put the r, u part 
into diagonal form 



ds^ = 
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-dr' 
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where t = f + f{u) with 



f{u) = 
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(10) 



which is akin to the trailing string solution. The fluc- 
tuations for u > Ui, are causally disconnected from the 
region u < u*, so that m* has precisely the form of a 
black-hole horizon. Euclidean regularity of this horizon 
imposes the temperature 
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(11) 



Restoring units |40| and assigning charge e = 1/vA to 
the fundamental charge carriers, as in [HI EH]) this may 
be expressed in terms of universal field theory quanti- 
ties as in Eq.Q. This effective temperature interpo- 
lates precisely between the two limits T* ^ ^/E and 
T^, ^ T which emerged from earlier calculations near 
the superconductor-insulator transition. Furthermore, 
the spatial fluctuations a = (a", a-*-) decouple from a* 
so that we are free to consider them separately. The 
fluctuations of interest satisfy the bulk equation 
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(12) 



where we have changed to the dimensionless variable z = 
m/u* and f{z) = 1 — z"**. The quantity Id := a;/(7rT*) is 
a dimensionless measure of frequency. We can now study 
the noise of the current j, dual to the bulk gauge field 
a, by following the prescription of [25] to calculate the 
Schwinger-Keldysh two-time correlators of the boundary 
field theory. The current fluctuations satisfy a Langevin 
equation of the form 



4) 

dt 



^{t,t')mdt' = m , 



(13) 



where we can extract the memory kernel j{t,t') — 
Gji{t,t') and noise ^(i) using holography, in analogy to 
holographic Brownian motion of heavy quarks [301 131| . 

By carefully studying the analyticity properties of pos- 
itive frequency modes on the Kruskal plane as we cross 
the horizon [53], we find the noise correlation 

(r(a.)e(-c.)) = G%Ju) = -(1 + 2n,)ImG'^;(u;) , (14) 

where G^,^{uj) is the symmetric/Keldysh correlator of 
the current ji in the Schwinger-Keldysh formalism and 
G^(a;) is the retarded correlator. We have the effective 
thermal factor n^, = l/(exp(w/r^,) — 1). This gives us a 
current noise 
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(15) 



where the final equality follows from a scaling analysis of 
Eq. ( 12 1. In fact, we can solve equation ( 12 1 exactly. The 



solution with ingoing boundary conditions at the horizon 
for both components of a is 



a.,{z) 



1 



z + 1 



exp 



ZtD 



arctan z 



(16) 



for constants q. 
Green function 



From this we deduce the exact retarded 
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(17) 



where cr(2-i-i) is the non-equilibrium conductivity deter- 
mined above. In turn this gives the current noise power 



Sj — 4ct(2+i)T* 



(18) 



The noise is related to the conductivity by an 
equilibrium-like fluctuation-dissipation or Einstein rela- 
tion at the effective temperature; an exact interpolation 
between the previously known results of W- 

Discussion: Hawking radiation is dual to quantum 
critical current noise under the holographic mapping. We 
have used this mapping in the probe brane limit to study 
non-linear current noise. The crucial feature of this anal- 
ysis is the effect of the electric field on the metric in- 
duced on the probe brane; this metric is modified such 
that a new horizon appears with which we may associate 
a temperature T*[4T]. This horizon describes the com- 
bined effects of thermal fluctuations from the equilibrium 
bath and pair production by the Schwinger mechanism. 
By considering the propagation of fluctuations from this 



horizon to the boundary, we have shown that the fluctu- 
ations on the boundary - and hence of the gauge theory 
- obey an equiUbrium-hke fluctuation-dissipation or Ein- 
stein relation at the effective temperature T^ . 

In general, an out-of-equilibrium steady state is not 
expected to obey an equilibrium-like fluctuation dissipa- 
tion relation. The present holographic analysis, suggests 
that this may nevertheless be the case in quantum crit- 
ical systems with a holographic dual. This rather ele- 
gant result emerges as a consequence of some very gen- 
eral properties of gravitational metrics. It would be very 
informative to re-cast it as a principle directly applicable 
to the condensed matter system. This principle is subtly 
hidden; it does not correspond to a thermal distribu- 
tion of the Klein- Gordon normal modes, since this would 
not carry a current. The distribution of fluctuations is 
a combined feature of the non-equilibrium distribution 
and the scattering integral. Nevertheless, thermal distri- 
butions of fluctuations in the presence of profoundly non- 
equilibrium distributions have been seen before in quan- 
tum critical systems in the low frequency limit j32[ I33j , 
even when they do not posses a known holographic dual. 
One might then speculate that if this principle can be 
identified, it will have broader application. 

It would seem that for a quantum critical system with 
a gravity dual, we may always consider non-linear re- 
sponse in the probe brane limit. Presumably, the driving 



field will always modify the induced metric in such a way 
that a new horizon is formed with which we may asso- 
ciate an effective temperature. As we have taken pains 
to demonstrate here, a fiuctuation-dissipation relation at 
this effective temperature will result. What are the lim- 
itations of this picture? We have used arguments from 
condensed matter to motivate the existence of a non- 
equilibrium steady state in the correct limit of thermal 
conductivity and system size. These conditions are not 
easily met in practice and certainly constrain the appli- 
cability of our results. Nevertheless, where they are met, 
it is unlikely that a gauge-gravity dual that applies in 
equilibrium would be invalidated by simply tuning up 
the strength of the probe field. One way to negotiate 
the constraints of system size and thermal conductivity, 
and an interesting avenue for further study, might be to 
study the approach to the steady state distribution in a 
quench after the sudden change of electric field similar to 
the mass and coupling quenches of j34j . 
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